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I. INTRODUCTION 



Neutrino oscillation physics has definitively entered the era of precision measurements 
of the fundamental neutrino parameters. Recent experiments such as Super-Kamiokande, 
SNO, KamLAND, K2K, and MINOS Q, 13, 3, 4, J 6|, have improved our knowledge of 
the neutrino mass squared differences {i.e., Amli and Arn^i) and some of the the leptonic 
mixing parameters {i.e., 612, ^23)- Whether 6*13 is different from zero and CP-violation is 
present in the leptonic sector of the Standard Model are questions which will be addressed 
by forthcoming experiments 

While there exist clear evidence that the V — A structure of the weak interactions of the 
Standard Model correctly describes neutrino interactions with matter, there is still a possibil- 
ity that some next-to-leading order mechanism affects the processes of neutrino production 
and detection. In general, this sort of Physics beyond the Standard Model is described by 
a set of higher- dimensional non-renormalizable operators suppressed by some high energy 
scale (see, e.q., Refs. and references therein). The precision measurements in the neu- 
trino sector then open up the possibility to investigate such non-standard interactions (NSI) 
at a quite accurate level. 

Neutrino oscillations and NSI in terrestrial neutrino experiments have been studied ex- 



tensive^ 
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y in the literature, using the neutrino factory project jiol . 
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20( 1 and other different neutrino facilities {e.g., conventional neutrino beams, super- 
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271] to assess the impact of the NSI in neutrino 



beams and /3-beams) 
physics. 

Here we adopt a different point of view, we investigate the NSI in the neutrino sector 
using very high energy neutrinos from astrophysical sources and the capability of neutrino 



telescopes to measure their fluxes on Earth 
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32|. We rely on the simplifled 



assumption that the new physics effects arise in the production and detection processes 
but do not affect the neutrino propagation. We also assume that the lepton mixing matrix 
is the standard unitary matrix describing the couplings of the charged-lepton-neutrino-l^ 
vertices. Moreover, for the sake of illustration, we prefer to study the possible NSI signals 
for the three different sources (pion, muon damped and neutron sources) separately, as the 
effects of new physics are quite different in these cases. 

This work is organized as follows. In Sec. [TTl we will present analytic considerations for 
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the NSI we study in the paper. In particular, we will first describe the detector effects in 
a unified way, since these are independent of the assumed neutrino source. We will then 
address the question of non-standard physics in the production processes. In Sec. Illlt we 
will describe the statistical approach we use to study the sensitivity of neutrino telescopes 
to new physics effects and also present our numerical results. Finally, in Sec. IIVI we will 
present a summary of the work as well as our conclusions. 



II. NSI AT THE SOURCE AND DETECTOR 



We consider NSI through effective four-fermion operators of the form 



B-NSI 



-2^PlGp cos(^h')£^J[m7''^lc?][47m^l^^/3] + h.c. 



(1) 



for charged-current NSI with baryons, where is the Weinberg angle, and 



-NSI 



-2 V2G^£^f %^^Pds,\ [u.i^Pl^s] 



(2) 



for purely leptonic NSI (see 33|, [SJ, [ssl, [36 1 for a discussion on the limits on While 
the addition of the hermitian conjugate in the case of NSI with baryons leads to operators 
of the form [(i7^PLM][%7^PL£a], the requirement of a hermitian Lagrangian for the purely 
leptonic NSI implies = ef"*. In what follows, we will assume that the coherence among 
the neutrino mass eigenstates that arrive at the Earth has been lost. Thus, we do not need 
to take any interference terms into account and may treat the source and detector processes 
separately. 



A. Detector effects 



Since the detector processes are essentially the same regardless of the source, we will start 
by considering these. In principle, the detection of astrophysical neutrinos is done through 
charged-current reactions of the form v + X Y + 1^ and then identifying the neutrino 
flavor through identification of the outgoing charged lepton. Assuming that the neutrino 
arrives at the Earth in the mass eigenstate i/j, the matrix element involved in computing the 
reaction rate is given by 

M,^ = M^[{l + e^'')UU, (3) 
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where Aio is the matrix element for the corresponding reaction if the incoming neutrino 
would have been of flavor a and there were no NSI, 1 is the 3 X 3 unit matrix, e"^ = (s^J), 
and U is the leptonic mixing matrix. Thus, if the incoming flux is (pi, then the measured 

flux of Ua is given by 

0, = 0i|[(l + £"')t/]a^|'. (4) 

Naturally, the composition of the neutrino flux arriving at the Earth is not purely Ui. 
Thus, in order to compute the actual flavor flux, we need to sum over all mass eigenstates 
and arrive at 

0a = ^0i|[(l+£"')t/]a^|'- (5) 
i 

The actual fluxes (pi are dependent on the source type and may also contain effects from the 
same NSI as those affecting the detector processes. 

B. Hadronic sources 

Let us flrst assume that the astrophysical neutrino source creates neutrinos through 
hadron decays only. Out of the usually considered scenarios, this category includes the 
muon-damped pion sources (giving initial flavor ratios of 0e ^ 0// ^ 0r = : 1 : in the 
neutrino fluxes) as well as neutron-like sources (1:0:0). These sources produce neutrinos 
through processes where a meson or baryon decays into a charged lepton of a given flavor^ , a 
neutrino and possibly other products. When a neutrino is produced together with a charged 
lepton of flavor /3 in a weak decay of a meson or a baryon, the probability of producing it 
in the neutrino mass eigenstate Ui is given by 

Pi = \Uf,i\' (6) 

in the Standard Model. With the addition of NSI, there is also the possibility to produce 
the neutrino in a different flavor state and the corresponding probability changes according 
to 

Pi^\[{l + s'''')UW. (7) 

^ In the case of pion decay, there is a small contamination of decays into electrons which is suppressed by 
ml/ml. 
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In general, the full expression for Pi also contains a normalization factor, since 1 + e'"'^ may 
not be unitary. However, this normalization factor is the same for all Pi and will be removed 
once we consider neutrino flux ratios. Since the meson/baryon decay gives the only neutrino 
contribution in this scenario, the Vi flux will be given by 

02 = (l)oPi, (8) 

where 0o is the total initial neutrino flux. By insertion into Eq. ([5]), the measured flux of 
at the detector is therefore given by 



oc 



(9) 



C. Non muon-damped pion sources 

We now consider the situation where the source is producing neutrinos both through an 
initial decay of a charged pion and through the subsequent decay of the resulting charged 
muon. The complete decay chain is then 

i , (10) 

as well as the corresponding CP-conjugate reaction for vr". With only Standard Model 
interactions, a = 7 = /i and f3 = e; the produced flavor ratio is then approximately 1:2:0. 

Mathematically, the decay of the pion can be described in the same way as in the previous 
section, with the exception that we now need to know the normalization factor for Pi, since 
we will add the contributions from the pion and muon decays and need to be consistent while 
doing so. The normalization factor A*"^ is simply given by the fact that the total probability 
of the pion decay should be equal to one, i.e., 

E^^= i^El[(l + ^"')f^Wl' = l =^ N, = [{1 + e-'){l + e-'%, (11) 

i ^ i 

(essentially, this is the factor by which the decay rate of the pion would change due to the 
NSI). 

The muon decay is a little more involved, since the final state involves two neutrinos, 
where one of these is not observed. In the literature, this problem is usually solved by 



considering only NSI of the form [e'y'^PLlAi^alpPL^e]- With this simphfication, the situation 
is completely analogous to the case of the pion decay, simply because we know the flavor of 
the outgoing V (in the case of /^+-decay). With general NSI however, we need to consider 
the full matrix element for the decay — > e^UiUj, which is given by 

Mij = Mo[U\j^' + enU]ji (12) 

where A^o is the matrix element when neutrinos do not mix and only Standard Model 
interactions are considered, J'^^ — Sea^ixp, and e^^^ — (s^^) is a matrix containing the 
strengths of the NSI. In order to obtain the probability for a neutrino from this type of 
decay to be in the mass eigenstate Ui, we need to consider the fact that we do not measure 
the anti-neutrino from the same decay. Thus, the probability will be given by an incoherent 
sum over the antineutrino mass eigenstates as 

Pt^ cx J2 °^ [U\J^' + enHJ''" + enU]u. (13) 

j 

Again, this needs to be normalized and the normalization factor is given by 

N^"^ = Tr[( J'^^ + £'^^)t( J/^^ + £^^)] (14) 
The corresponding argumentation for the outgoing antineutrino results in 

Pf^ oc J2 « [U\J^' + s^'){J^' + en^U]jj (15) 

i 

with the same normalization constant iV^^ = A^'^^. Repeating the same derivation for /x^- 
decay, we arrive at P/* — Pj^^ and Pj* — P^^ . Thus, assuming equal numbers of positive 
and negative pions decaying, the flux of neutrino mass eigenstate (or antineutrino mass 
eigenstate Pj) is given by 

0, = (/-o(P, + /^'^%P,'^^), (16) 
where 0o is the initial flux of Ve in the case of Standard Model interactions only. 



III. EXPERIMENTAL IMPLICATIONS 



For experimental reasons, it is customary to consider the flavor flux ratios 

Re, = ^. R,r^^, and (17) 



(note that these are not independent). With the flavor fluxes computed as in Sec. [TTl we want 
to know how NSI at the source and detector can affect the results of these measurements. 
In order to do this, we will examine how well a measured value of any of the Rs can be 
accommodated when considering Standard Model interactions only, as well as when allowing 
for NSI. To quantify this, we introduce the function 

x^w (18) 

Here, i?^^^ is the measured value of a given ratio, T^^[x) is the theoretically expected value 
computed from the set x of input parameters, ctr^ is the experimental uncertainty in 
i is an index running over the flux ratios considered, a^^ are the external uncertainties for 
the input parameters, and j is an index running over the input parameters. In the case 
of Standard Model interactions only, x contains the neutrino mixing parameters, while in 
the NSI scenario it also contains all of the different es. How well a particular measured 
combination of W^^^s can be accommodated within a model is then quantified by 

xLin = min[x^(x)], (19) 

which is distributed with n degrees of freedom, where n is the number of independent 
flavor flux ratios considered. 

In the following numerical com put ations, we adopt the standard parametrization of the 



371 . l38l . l39l . l4d . l4l| ) , assuming also a tri-bimaximal struc- 



neutrino mixing matrix (see Refs. 
ture 



421 for it. We assume the 3a errors on the standard mixing parameters to be of the 

n 

same order as the errors in Ref. m. In Tab.m we give the corresponding values for these 
parameters. 

A. Muon-damped pion sources 

In order to appreciate the effects of the new physics on the flux ratios of Eq. ( ITTl) . it is 
useful to recall the Standard Model expectations for these quantities. Putting e^'^ = and 
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Parameter 


Best fit 


3(7 error 




35.3° 


5° 


^13 





12.5° 


^23 


45° 


10° 


6 


Free ([0,27r]) 



TABLE I: Summary of the best fit values and 3a errors used in our numerical computation for the 
standard neutrino oscillation parameters. 

f3 = fi in Eq. we get the following expressions for the flavor fluxes: 

^ ^ 1 2 1 1 2 

0e = 00 > \U,,i\ lUpj] 



i 
i 



4 

fll I ; 



(20) 



In the limit of exact tri-bimaximal mixing 
we obtain: 



^13 



0, sin^^is = 1/3, and sin^ ^23 = 1/2 |42j), 



4 
7' 



7 

and R = —. 



(21) 



Before concluding that large deviations from these values are signals of new physics, we 
should carefully take the role of the uncertainties of the standard parameters into account 



44( 1 . In order to do that, we compute the distributions of the three ratios, extracting 
neutrino mixing angles according to Gaussian distributions with central values and la error 
as deducible from Tab. [B The result of this procedure is shown in Fig. [H It can be clearly 
seen that the larger spread is obtained for i?e^, due to the uncertainties on the mixing angles 
6'i2; ^23 and on the product cos 5^13 (see also Eq. (4) in Ref. 45|]), which are all of the same 
order of magnitude. On the other hand, the fact that most decays will result in neutrinos of 
the mass eigenstate with the largest content protects from becoming much smaller 
than one, an effect that can also be seen in the ratio R, which cannot be much smaller than 
1/2. The pile-up of the distributions close to these values are the results of 6^23 being close 
to maximal. It is then possible that new physics effects would be much more visible in the 
ratios containing r neutrinos and, for this reason, we will investigate the effects of the e 
parameters in the i2^T--i?-plane. 
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FIG. 1: Statistical distribution (in arbitrary units) of the flux ratios Re^ (dashed line), R^r (solid 
line) and R (long-dashed line), as computed in the Standard Model for a muon-damped pion source. 



The correlation between i?^^ and R can be analyzed using the definition in Eq. f|T8l) . 
A crucial point here is the precision we expect in the experimental measurement of the 
fiux ratios, a^-; for the sake of illustration, we will assume that a^- = O.lRf^^ for any 
ratio (the effect of changing this value is discussed in Sec. IIII El) , which means that they 
are measured with a 10 % error. In addition, we assume an external error of 0.1 for the 
e parameters. The reason for choosing such a large error is that we are mainly interested 
in the qualitative impact of the parameters rather than making precise predictions. The 
results of the minimization of the over all model parameters is shown in Fig. [2J We 
note that the extension of the isocontours produced by new physics is predominantly in the 
direction of large ratios; in particular, R^^r and R can be as large as twice their tri-bimaximal 
mixing values, thus an experimental signal in this direction could be ascribed to the effects 
of non-standard interaction type described in this paper. The extension of the contours to 
small values of the flavor flux ratios {e.g., R^r ~ 0.5, to be compared with Fig. [1]) even in 
the standard case can be attributed to the flnite resolution assumed for the measurements. 

It is interesting to understand which of the new parameters that cause the largest de- 
viations from the Standard Model prediction. This is illustrated in Fig. [3], where we show 
the isocontours of Xmm -R^r^-R-plane when allowing for only one non-zero e at a time. 
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FIG. 2: Isocontours of Xmm R^r^R-pldiT'^ the case of a muon damped pion source in 

the standard framework for neutrino oscillation and interaction (shaded regions) as well as for the 
framework with included NSI (black curves). The contours correspond to 90 %, 95 %, and 99 % 
C.L., respectively. The black cross corresponds to the prediction in the case of tribimaximal mixing 
in the standard framework. 

From this figure, it is clear that the most relevant contributions to the extension come from 
£^'lj_,£^f and e^f. Thus, these are the parameters to which neutrino telescopes would be 
most sensitive. Putting stringent bounds on these parameters would therefore be useful for 
making more robust predictions for the flavor flux ratios in the case of the muon-damped 
pion source. Note that none of the outer contours from varying individual e reach the 99 % 
C.L. contour from the full simulation. This implies that in order to accommodate these 
contours in the full simulation, some combination of different es is needed. 

B. Neutron-like sources 

In the case when the source of the astrophysical neutrino flux is neutron-like {i.e., es- 
sentially a beta decay), then the considerations corresponding to those made for the muon- 
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FIG. 3: Isocontours o/ Xmin R^T^R-p^O''ne in the case of a muon- damped pion source. The 

black curves correspond to the case with only Standard Model interactions, while the three inner 
shaded regions are produced allowing one e to vary with an external input error of 0.1 at la. The 
panels are placed such that the e varied has the same position in the e"""^ matrix as the plot has in 
the figure, i.e., the upper right plot corresponds to allowing e"^ to vary (we put the neutrino flavors 
in the order {ve^v^^Vr})- The contours correspond to 90 %, 95 %, and 99 % C.L., respectively. 
The outer shaded region is the 99 % C.L. contour from the full NSI simulation. The black crosses 
correspond to the prediction from tribimaximal mixing in the standard framework. 



damped pion sources result in the flavor flux ratios 

i?eM = | Rf.r = l, and i? = ^, (22) 

in the case of tribimaximal mixing. Figure H] shows the isocontours of Xmm the R^r^R- 
plane for a neutron-like source using the same procedure as that described in the previous 
section. In this case, the extension of the isocontours is mainly in the direction of increasing 
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FIG. 4: Isocontours o/Xmin Rfj_r-R-plane in the case of a neutron-like source in the standard 

framework for neutrino oscillation and interaction (shaded regions) as well as for the framework 
with included NSI (black curves). The contours correspond to 90 %, 95 %, and 99 % C.L., re- 
spectively. The black cross corresponds to the prediction in the case of tribimaximal mixing in the 
standard framework. 

Rfj,r, corresponding in an increase in the z/g flux on the expense of Ur compared to the 
standard setup. However, this extension more or less has the same shape as the extension in 
the standard case and the big errors for large values of the flavor flux ratios can be somewhat 
attributed to the increasing absolute errors on a/?.. 

In Fig. [5], we show the impacts of the individual e parameters for a neutron-like source. 
Unlike the case of the muon-damped pion source, with a neutron-like source we find that all 
of the 6 has at least some small impact on the flavor flux ratios, even if it is still relatively 
small. The largest difference is observed for the e"^ parameter, which (together with e"^) 
seems to be the source of the extension to larger i?^^ with R being essentially constant. Also 
the impact of e^^ is peculiar, even if it is not as large as that of e^J^. It would seem that this 
parameter fixes i?^,- while altering R, signifying a change in the flux of z/g compared to the 
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FIG. 5: Isocontours o/ Xmin Rfj^r ^R-P^me in the case of a neutron-like source. The black 

curves correspond to the case with only Standard Model interactions, while the three inner shaded 
regions are produced allowing one e to vary with an external input error of 0.1 at la. The panels 
are placed such that the e varied has the same position in the e^'^ matrix as the plot has in the 
figure, i.e., the upper right plot corresponds to allowing e^^ to vary (we put the neutrino flavors 
in the order {V(,,i'^,i't})- The contours correspond to 90 %, 95 %, and 99 % C.L., respectively. 
The outer shaded region is the 99 % C.L. contour from the full NSI simulation. The black crosses 
correspond to the prediction from tribimaximal mixing in the standard framework. 

other two species. 

C. Non muon-damped pion sources 



As discussed in Sec. [TTl when a pion source is not muon-damped, there will essentially be 
two different sources of neutrinos - the pion decay and the subsequent decay of the muon. 
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FIG. 6: Isocontours o/Xmin R^r^R-plo-ne in the case of a non muon- damped pion source in 

the standard framework for neutrino oscillation and interaction (shaded regions) as well as for the 
framework with included NSI (black curves). The contours correspond to 90 %, 95 %, and 99 % 
C.L., respectively. The black cross corresponds to the prediction in the case of trihimaximal mixing 
in the standard framework. 

Since the NSI involved in the two different processes do not depend upon the same es, the 
effective parameter space is increased quite dramatically, even compared to the previously 
considered NSI scenarios. For this type of source, the prediction of tribimaximal lepton 
mixing is given by 

= 1, Rf,r = 1, and = (23) 

simply due to the fact that the flavor fluxes are predicted to be equal. 

In Fig. El we show the results of the minimization of both for the standard case and 
when allowing for NSI (both in the pion and muon decays). For large values of R^r and R 
in the NSI case, this figure is not very smooth. The reason for this is that the is now a 
function of a very large number of parameters and the minimization procedure is not always 
able to find the actual minimum. In this scenario, we see that the extension of the contours 
due to NSI is mainly along the direction of constant R/ R^r-, corresponding to a fixed ratio 
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between the Vf. and fluxes while the relative Vy, flux varies (again, the smaller size of the 
regions for smaller Ri is mainly due to the differences in cxj?.). 

As in the previous scenarios, we present the dependence on individual es in Fig. [71 Here, 
the exception is that the parameter space now consists of an additional nine es, which we 
illustrate by two different plots for the dependencies on the e^'^ and 5^*^, respectively. 

The effects of the e^'^ parameters are very similar to what they were in the case of the 
muon-damped pion source, although a bit less pronounced. This is to be expected, since the 
processes of the muon-damped source also play a role in the case when the muons are not 
damped - although the effect is less pronounced as there are also other processes involved. 
However, the effects of the e^^ parameters are essentially negligible, with the only parameter 
actually showing any difference (although still a very small one) being e^^. Thus, it would 
seem that the cannot play a signiflcant role when it comes to astrophysical neutrino 
fluxes. 

D. Comparison of scenarios 

An important aspect in the study of astrophysical neutrino fluxes is what we could learn 
about the source from studying the flavor composition of the flux (z.e., what type of source 
that is providing the flux). Thus, in Fig. [HI we show the 99 % CL. contours for all of 
the source types discussed above. As can be seen in this flgure, the pion sources could 
be relatively well separated from the neutron-like sources in the standard scenario, while 
it could be harder to tell whether there is muon-damping or not in a pion source. This 
conclusion is not signiflcantly altered when taking NSI into account, although it will be a 
little harder to tell the sources apart if flavor ratios in the region expected from Standard 
Model interactions only are measured. However, since the NSI mainly extend the contours in 
different directions for different scenarios, large parts of the space of possible measurements 
allowed in the case of a neutron-like source are not allowed for pion sources and vice versa. 
Still, telling if there is muon-damping in a pion source remains a difficult endeavour. The 
contour extension due the NSI also has another intriguing property. Since the extensions 
are mainly not in the direction to where other scenarios are located, there is a fair chance 
that if NSI actually appear in the flavor flux ratios, their effect could not be misinterpreted 
as simply being due to a different neutrino source. 
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FIG. 7: Isocontours of Xmin R^T^R-plo,ne in the case of a non muon-damped pion source. 

The black curves correspond to the case with only Standard Model interactions, while the three 
inner shaded regions are produced allowing one e to vary with an external input error of 0.1 at la. 
The upper plot shows the dependence on the e"*^ and the lower plot the dependence on the e^'^ . The 
panels are placed such that the e varied has the same position in the matrices as the plot has in 
the figure, i.e., the upper right panel of the upper plot corresponds to allowing e^'^ to vary. The 
elements of the figure are the same as in the previous plots. 
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FIG. 8: The 99 % C.L. isocontours in the R^r^R-plane for the sources we have discussed. The 
solid lines correspond to the contours when allowing for general NSI, while the dashed contours 
correspond to the result in the case of Standard Model interactions only. 

E. Impact of experimental uncertainties 

As mentioned earlier, the experimental precision with which the neutrino flavor flux ratios 
can be measured is an important part of our analysis. In our discussions so far, we have 
assumed a relative error of 10 % and it is of importance to understand how this affects our 
results. For this purpose. Fig. [9] shows the full result for the muon-damped source both 
with our initial assumption of = O.li?^^'' as well as with a more optimistic assumption 
of cxr. = 0.05Rf^^. As can be seen in this figure, a better experimental resolution makes 
the NSI degeneracies more pronounced. This is to be expected, since better measurements 
allow for better opportunities to extract information of the underlying physics. Conversely, 
if the resolution is made worse, it will become harder to distinguish the NSI degeneracies 
from the standard scenario as deviations could simply be due to experimental uncertainties. 
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FIG. 9: The isocontours in the R^r^R-p^o-ne in the case of a muon-damped source assuming differ- 
ent experimental resolutions for the flavor flux ratios. The left panel assumes a^. = O.li?^^^ (^i-e., 
identical to Fig. while the right panel assumes aji- = 0.05i?^^^. The elements of the plots are 
the same as in Fig. [H 

IV. SUMMARY AND CONCLUSIONS 

In this paper we analyzed the impact of new physics effects in the production and de- 
tection of astrophysical neutrinos. Detection processes were considered independent of the 
source producing neutrinos and the possible new physics effects were parametrized in terms 
of the matrix e'^'p which alters the charged- current reaction u + X — > Y + ia involving 
leptons and hadrons. The same matrix also intervenes in the production of neutrinos from 
muon-damped pion sources (in which neutrinos coming from subsequent muon decays do not 
contribute to the measurable neutrino flux on Earth) and from neutron-like sources. If muon 
decays also play a role in the production of neutrinos (as in the case of the "standard" pion 
sources), then a new set of parameters e^^^ should be introduced to take the purely leptonic 
process into account. Since the fluxes of astrophysical neutrinos can be measured at neutrino 
telescopes, we have considered the experimentally accessible flux ratios R^^r = 'P^/'Pt and 
R = 4'fj,/{(f)e + 4>t), in order to investigate the effects of the new parameters. This was done 
through a analysis to see how well the neutrino flux ratios at neutrino telescopes can 
be accomodated in different models, assuming tri-bimaximal mixing as an external central 
value for the standard neutrino oscillation parameters. Making this analysis in the Stan- 
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dard Model first and then including the NSI gave us insight of how new physics can affect 
the possible ranges of these observables. Our results can be summarized in the simplified 
scenario, where only one source at a time is responsible for the neutrino flux on Earth. In 
particular, for muon-damped sources, we found that with NSI, the measured flux ratios can 
be significantly larger than their tri-bimaximal mixing values, an effect mainly due to eJJ^, 
and £"^. Also for neutron- like sources the extension of the contour plots obtained including 
NSI effects are in the direction of increasing R^r, with R being essentially constant (thus a 
significant change in the flux is at work). Moreover, all the epsilon parameters contribute, 
in some extent, to the shift of the flavor flux ratios. The analysis of the non muon-damped 
sources is more complicated. This is due to the two independent classes of new parame- 
ters entering in the game. However, while both R and R^r can be sizably different from 
the Standard Model predictions, we found that the leptonic e'^^ have negligible impacts on 
these results. In principle it is possible, in both Standard Model and new physics scenarios, 
to learn something about the source, since the observational parameter space covered by 
the flux ratios from one particular source does not overlap significantly with that of another 
source. 

It should be noted that all of our considerations depend on the assumed uncertainties 
on the extraction of the flux ratios as well as on how well the standard neutrino oscillation 
parameters have been determined. In fact, although smaller ctr. makes the NSI effects 
more pronounced, larger (TR^ will completely oversahdow any significant new physics effect. 
In addition, we have considered relatively large values of the NSI paramters. It should 
therefore also be noted that the effect of small NSI will become increasingly difficult to 
observe unless the flavor flux ratios and standard neutrino oscillation paramters can be very 
accurately determined. 
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